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We study equilibrium states of quantum spin systems with nonadditive long-range interactions 
by adopting an appropriate scaling of the interaction strength, i.e., the so called Kac prescription. 
In classical spin systems, it is known that the equilibrium free energy is obtained by minimizing the 
free energy functional over the coarse-grained magnetization. Here we show that it is also true for 
quantum spin systems. From this observation, it is found that when the canonical ensemble and the 
microcanonical ensemble are not equivalent in some parameter region, it is not necessarily justified 
to replace the actual long-range interaction by the infinite-range interaction (Curie- Weiss type in- 
teraction). On the other hand, in the parameter region where the two ensembles are equivalent, this 
replacement is always justified. We examine the Heisenberg XXZ model as an illustrative example 
and discuss the relation to experiments. 
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I. INTRODUCTION 

Long-range interacting systems have attracted much 
attention because of a wide variety of physical systems 
and peculiarities of its dynamics and thermodynam- 
ics [1, 2]. Recent experimental achievement allows us to 
access artificial quantum lattice systems where the cou- 
pling parameters can be controlled. Dipolar gases in op- 
tical traps are suggested to be possible experimental re- 
alization of long-range interacting lattice systems. More- 
over, it has been suggested [3] that we can realize artifi- 
cial gravitational systems in laboratory by irradiating off- 
resonant laser beam onto atoms in a Bose-Einstein con- 
densate. Theoretically, long-range interactions exhibit 
rich equilibrium and non-equilibrium behavior like the 
ensemble inequivalence, negative specific heat, and the 
absence of thermalization. These properties are atypical 
for short-range interacting systems. Therefore, studying 
long-range interacting systems will lead us to deepen our 
fundamental understanding of statistical mechanics. 

Those peculiar features of long-range interacting sys- 
tems have been understood by using mean-field mod- 
els [1]. In mean-field models, the interaction between par- 
ticles or spins is assumed to be independent of distance. 
Since exact results are available in mean-field models, 
they are important theoretically although they appear 
to be too simplified. Thus exploring the condition on the 
realization of predictions by mean-field models should be 
actively proceeded. 

Kastner [4, 5] examined quantum spin systems with 
an infinite-range ferromagnetic interaction. He demon- 
strated the ensemble inequivalence of these systems. 
In this paper, more general long-range interactions are 
treated; ferromagnetic interactions whose range is com- 
parable with the system size. An important class of these 
interactions is the power- law potential 4>(r) ~ —l/r a , 
where r is the spatial distance between spins and a sat- 
isfies < a < d (d is the spatial dimension). It has been 
believed from the analysis of exact partition functions or 
the numerical calculations of some models that thermo- 



dynamic behavior does not depend on a as long as it is 
less than d [6-9] . As a = corresponds to the mean-field 
model, it implies that the results of the mean-field model 
might be exact for general long-range interacting systems 
as long as <f>(r) ~ l/r a , and < a < d. 

In the previous papers [10-12], it was shown that for 
classical spin systems, the free energy does not depend 
on a in the canonical ensemble with non-conserved order 
parameters, but a becomes relevant in some parameter 
region when there is a constraint of the value of the en- 
ergy or the magnetization. A sufficient and a necessary 
condition for a to be relevant are derived in the previ- 
ous works. In this paper, we show that it is also true in 
quantum spin systems. 

An immediate conclusion is that if the inequivalence 
of the canonical and the microcanonical ensembles pre- 
dicted by the analysis of the mean-field model can be 
observed in experiments, the violation of exactness of 
the mean-field theory will be also observed in the micro- 
canonical ensemble. Inhomogeneous spin configurations 
in equilibrium are its sign. Although it is known that 
the mean-field models reproduce many properties of gen- 
eral long-range interacting systems even quantitatively, 
the infinite-range (or Curie- Weiss type) interaction does 
not become an adequate idealization of actual long-range 
interactions in such a situation. 

This paper is organized as follows. In Sec. II, we in- 
troduce the model. In Sec. Ill, an important notion of 
coarse-graining is explained. In Sec. IV, we express the 
free energy only by classical variables. In Sec. V, the 
variational expression of the free energy is derived. In 
Sec. VI, we mention some results deduced from the pre- 
vious sections, namely that the free energy is independent 
of the precise form of the interaction potential as long as 
there is no constraint such as fixed value of the energy or 
the magnetization. On the other hand, the free energy 
under the constraint of the fixed value of the magneti- 
zation and the microcanonical entropy may depend on 
the interaction potential. This result is the same as in 
classical spin systems with long-range interactions. In 
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Sec. VII, we investigate the spin-1/2 XXZ model in the 
canonical and the microcanonical ensemble as an illus- 
trative example. In Sec. VIII, we summarize our results. 

II. MODEL 

We study systems expressed by the following Hamilto- 
nian on d-dimensional regular lattice: 

n r n 

tf = -^E^(^) e - r 5> « 

ij a—x.y^z i 

We focus on the ferromagnetic interaction, A a > and 
(/)(■) > 0. The parameter 7 corresponds to the inverse of 
the interaction range and we take the limit of 7 — > ap- 
propriately (see below). </>(•) is the interaction potential 
and we will give the condition of <f> later. s*j is the spin 
operator at site i whose length s is independent of par- 
ticular site i and s e {1/2,1,3/2,...}. Normalizations 
1 /s 2 for the interaction strength and 1 / s for the magnetic 
field h are not essential. With these normalizations, the 
system is reduced to a classical continuous spin system in 
the limit of s — > 00. The lattice interval is set to be unity 
and L denotes the length of the system (the number of 
spins is L d = N). The position of site i is denoted by 
Ti 6 Z d n [1, L] d . The distance between two sites i and j 
is denoted by r^. 

We take the two limits; L — > 00 and 7 — > 0. In the 
present paper, we consider either of the following limiting 
procedures: 

(i) nonadditive limit: the condition for (j>(r) is that 
4>(r) > and J r£ [ 7£ ]d (f>{r)d d r = 1. We take the 
limit 7 — > 0, L — > 00 with jL fixed. Note that 
it includes power- law potentials </>(r) ~ l/r a , < 
a < d. 

(ii) van der Waals limit: the potential (j)(r) satisfies 
4>{r) > and J <fi(r)d d r < +00. We take the van 
der Waals limit, that is 7 — > after L — > 00. 

The case (i) corresponds to nonadditive interactions, 
namely the interaction range is comparable with the sys- 
tem size. On the other hand, in case (ii), the interaction 
range is much longer than the lattice interval but much 
shorter than the system size L. Hence, two limiting pro- 
cedures treat different situations. 

Lebowitz and Penrose [13] proved that in the van der 
Waals limit, the equation of state of the system is equal 
to the mean-field equation of state (van der Waals like 
equation) with the so called Maxwell's equal area rule (or 
referred to be as the Maxwell construction). In the case 
of (ii), therefore, the free energy is independent of the 
precise form of the interaction potential </>(•). Extension 
to quantum many particle systems was done by Lieb [14] . 
The perturbative analysis around 7 = was suggested to 
understand the properties of a system with long but finite 



range interaction (7 > 0) [15]. Even in the critical region, 
some exact results were obtained along this line [16]. 

On the other hand, the case of (i) covers much longer 
interactions than the van der Waals limit. Similarly to 
the case of (ii), it has been indicated that the free en- 
ergy does not depend on the interaction potential, and 
it is equivalent to the mean-field free energy without the 
Maxwell construction in the nonadditive limit. It is called 
the exactness of the mean- field theory [6]. Therefore, 
we can regard the exactness of the mean-field theory in 
the nonadditive limit is a generalization of the Lebowitz- 
Penrose theorem in the van der Waals limit. Actually, 
it is correct in classical spin systems unless there are 
some constraints such as fixed energy or fixed magne- 
tization [10-12], as mentioned in the Introduction. 

III. COARSE GRAINING 

In this section, we introduce a kind of coarse grain- 
ing and show that the equilibrium free energy is exactly 
given by the free energy of the coarse-grained model. We 
divide the system into many cells which are large enough 
to contain many spins, but small compared to the in- 
teraction range. The length of each cell is denoted by 
I and cells are labeled by C p , p — 1,2,..., (L/l) d . The 
resultant spin of a cell C p is denoted by S p = X)iec 

Eigenvalues of (S p ) 2 and are denoted by J P (J P + 1) 
and M p , respectively (we set h to be unity). In each cell 
C p , there are l d spins, so they are properly labeled by 
s p (i), s p (2), . . . , s p (i<iy We define n-partial resultant spin 
as 

n 

^ n) =E^W. ( 2 ) 
i=i 

and eigenvalues of (Sp^) 2 arc denoted by J p n \j p n ^ + 1). 
Then the full Hilbert space is spanned by the ten- 

sor product of eigenstates \J p ,M p ,T p ') of each cell, 

® p L Ji )d \Jp,M p ,T p ld) ). Here 

r(i d ) = jo 7« /( 2 ) /C") - / \ 

is the "trajectory" of J p n ^ , n — 1, 2, . . . , l d . It is explained 
later in more detail. 

When the interaction is long-range, it is expected that 
the Hamiltonian can be expressed only by the resultant 
spin operators {S p } approximately. Moreover, typically 
J p is considered to be very large when I 3> 1, so it is 
naturally expected that all the S p , p = 1, 2, . . . , (L/l) d , 
can be regarded as classical spins whose length is J p . 

We shall construct the above argument rigorously. The 
coarse-grained Hamiltonian is introduced by 

(L/l) d - (L/l) d 

h = -^J2 i d <t>p q E s P sa q -- s - E 4 ( 3 ) 

p.q a—x,y,z p 
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where 



2d E E 0(7T«). 



(4) 



ieCp jec q 



The above coarse-grained Hamiltonian is expressed only 
by the resultant spin operators of each cell. The original 
Hamiltonian is well approximated by this coarse-grained 
one. In order to see it, consider the operator 



We define the norm of the operator A by 
\\A\\= sup |<VW)|. 

The norm of g then is given by 
L d 



(5) 



(6) 



2 7 <V 



(L/l) d 

E E E w™) 



E 



4.E EEi«- 



'pgl 



EA s a s a 



n 2 



Here ||£ a A a a?s? || /s 2 < A n 
fore 



(7) 

max a=a;i j /i2 A a , there- 



< 



E E E w™) 



Sel- 
ls) 



p,q iec p jec q 
In Ref. [11], it is shown that 

lim lim g c \ = 

l— too L— >oo 

for the nonadditivc limit and 



lim lim lim g c \ = 

for the van der Waals limit. Therefore, the original 
Hamiltonian is well approximated by the coarse-grained 
one. For convenience, we use the same notation "Lim" 
for the nonadditive limit and the van der Waals limit. 
The above property then implies 



LimlMI = 0. 



(9) 



By using this property, we shall show that the free 
energy per spin is exactly the same as that given by the 
coarse-grained Hamiltonian. The free energy is defined 

by 



(10) 



Here f3 = 1/T is the inverse temperature. The coarse- 
grained free energy is given by 



1 



mh) = -j^lnTre-^. 



(11) 



From now on, we show that Lim|/ — /| =0. We use 
the Bogoliubov-Peierls inequality, 



(12) 



where {\i)} is arbitrary orthonormal set. From H = H + 
L d g and Eq. (12), we obtain 



<i|fll<> 



(13) 



By choosing the eigenstates of H as {\i}}, we have 
Tre-? 11 > ^-/^-/^(ilsl*) > e^^Tre-^, 

(14) 

where {Ei} are eigenvalues of H. 

On the other hand, the upper bound of the parti- 
tion function Tre~^ H is obtained by using the Golden- 
Thompson inequality [17, 18], 

Tre"^ = Tre-P 6 -^'* < Tre^V^Y (15) 
We have 

Tre-^e"^ = E e-^(i\e~^\i) 

i 

<^e-^||e-^|| 



e /3£ d ||3ll Tre -^_ 



From Eq. (13) and Eq. (16), we obtain 

Tre-^ g 



e -PL d \\g\\ < Tre m < e PL d \\g\\ 



and thus 



l/-/l<IMI- 



(16) 



(17) 



(18) 



FromEq. (9), Lim|/— f\ = 0. Hence we can safely replace 
H by H to compute the free energy. 

IV. CLASSICAL BOUNDS 

The coarse-grained Hamiltonian depends only on the 
resultant spins {S p }. Since there are l d spins inside a cell 
and finally we take the limit of I — > oo, it is expected that 
the typical length of resultant spins is huge. We are then 
able to replace the quantum mechanical operator S p by 
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the classical vector. In this section, we develop a theory 
based on the above argument rigorously. 

We consider the Hilbert space of a cell C p . The 
spins in this cell is labeled by s p (j), i — 1,2, . . . ,l d . 
As is mentioned in Sec. Ill, the Hilbert space of the 
cell C p is spanned by the complete orthonormal basis 

{\j p ,M p ,rp)}. 

We notice that when each spin Si is a spin-s operator, 

e {4 k -^ + s, 4 k -^ +5-1, 4 k - 1] +8-2, 
....J^-aJniO.oo). 



(19) 



The "trajectory" with length l d is defined by 



where {Jp^} satisfies Eq. (19). We will omit (0 -> J p ) 
when there is no confusion. 

The set of trajectories of all T p ' (0 — > J p ) with a fixed 
J p is denoted by T[d(J p ). Moreover, we define the weight 
function W{J p ) by 



W(J P ) = [the number of elements of Ti<t(J p )] . (20) 

Explicit calculation of the weight function is given in Ap- 
pendix A. 

The identity operator in the Hilbert space of a cell C p 
is expressed by 



i P = E E E \J P ,M p ,rp) {Jp ,M p ,rp\. 

J p M p =-J Pr v* )f . Ttd(Jr) 

(21) 

Because the coarse-grained Hamiltonian depends only on 
the redundant spins S p , the partition function of the 
coarse-grained Hamiltonian is written as 



(L/l) d 

Tre-^= E E 



P =i 

(L/i) d 

® 

P =i 



(L/Z) d 



(j p ,M p ,rf)| 



-/3if 



|J P ,M p ,rf)) 



P =i 



(L/0 d 

n \y, e *w 

p=l l J v Mv = -.J„ 





"(L/0 d 




"(L/0 d 


X 






kp,M p 




p=i 




P =i 



(L/0 d 

- n (e^)|v^ 
p=i i j„ 



(22) 



Here |J p ,M p ) is an arbitrary state vector such that 

\j p ,m p ) e {|J p ,M p ,rf )),rf ) e 71, (J P )}. 

For instance, we choose |J p ,M p ) as 
1 



|J P ,M P ) 



E \J P ,M p ,T p l >>. (23) 



rj,' 'er ld (J P ) 



The symbol Tr{ Jp j = tr^tr,^ . . . tej {L/l)d means the trace 
over the subspace of fixed values of { J p }, 

tr Jp A = (J p ,M p \A\J p ,M p ). (24) 

M p = -J p 

Lieb [19] proved that for the general Hamiltonian such 
as Eq. (3), if (f> pp > 0, 

ZCl{{Jp}) -^^)- ZCl{{Jp + 1}l (25) 

where Z cl ({J p }) is the classical partition function with 
spins {J p }, namely, 



(L/iy 



x cxp 



47T 



-f3H c \{J p ,e p ,<P p }) . (26) 



Here, the classical Hamiltonian is defined by the same 
form of Eq. (3) but the spin operators S p are replaced by 
the classical vectors, 

S p — > J p (sin 9 P cos </> p , sin 9 P sin <j> p , cos 9 P ). 

Lieb proved the inequality (25) by using a spin coherent 
state representation. 

From Eqs. (22) and (25), we obtain the classical 
bounds 



z_< z <z. 



+ 1 



(27) 



where 



z - = E w a j p}) f m 
{jp} 



\{{2J P + 1) 



xcxp [-I3H C \{J P ,9 P ,^ P }) 



{Jp} 



x exp 



}J(2J p + 1) 



[-PH C \{ J p + l,0 P Ap})] 



Here W({J P }) is defined by W({J P }) = Y[ p W{J p ), and 
Df2 = dflidfl2 ■ ■ ■ dCl( L /i}d. 



■5 



V. VARIATIONAL EXPRESSION OF THE 
FREE ENERGY 

In this section, we evaluate the upper and the lower 
bounds of the free energy and derive its variational ex- 
pression. Let us define 

inf H c \{J p ,9 p ,<t> p }) = E*({J p }). (28) 

{Vp,<pp} 



We then find 

{J P } 



\L/l) d 

p 



-fiE*({Jp + l}) 



(29) 

Because £ {Jp} ./({ J P }) < {2sl d )W max {Jp} /({ J p }) 
for an arbitrary function f{{J p }), 



Z+ < (2sl d ) (L/l)d expl max 



pE*({J p + l}) 



{L/iy 



+ ln(2J p + 1) + In W({J P }) 



(30) 



Therefore, we obtain 



Lim/ > Lim inf — r 

{s P } L d 



H c \{j p + i,e p ,<j> P }) 



--glnW({J p }) 



(31) 



We show that we can replace J p + 1 by J p in Eq. (31). 
It is verified by the following evaluation: 



h c \{j p + 1, e p , <f> p }) - h c \{j p , e p , <j> p }) 

Jp + Jq + l 



2s 2 



(L/ir 

E ^ E 

pq a—x,y,z 



(32) 



where e p — S P /\S P \. As J p < sl d , we have 

h c \{j p + 1, e p , P }) - h c \{ j p , e p , P }) 



2sl d + 1 



(L/iy 



< 2g2 A max 2^ 7 < 



P9 



Here we note that 

(L/l) d 



E ^ = ^E7^(7n,) = ^, (34) 



(33) 



Thus we obtain the lower bound 



Lim/ > Lim inf — 7 

{5 P } ^ 



iF 1 ({J p ,# p ,0 p })--lnW({J p }) 



Next we estimate the upper bound. We define 



(36) 



(9* p , 4>* p ) = arginf H c \{J pi 9 P , (37) 

{0 P ,<Pp} 

9' p = \8 p - 9*1, and <$>' = mod((/> p - <f>*, 2w). In addition, 
the set w p (A) is defined by 



Wp (A) = {(0 p ,0 p ):tf p e[O,A]}. 



(38) 



For sufficiently small A, the volume of the set w p (A) 
satisfies 



1(6 



d% = 2tt(1 - cos A) > A 2 . (39) 



0p,0p)ew P (A) 
Moreover, we can show that 

sup H c \{J p ,9 p ,(b p }) -E*({J P }) < L d A max A, 
(6 P ,4>p)euj p (A) 

(40) 

whose proof is given in Appendix B. 
From Eqs. (39) and (40), 



xexp [~pH c \{J p ,9 p ,<p p }) 



x exp 



(47r)( L /') d 



-0 sup H c \{J pi 9 pi <j> p }) 
{(0p,0 p )ew p (A)} 



/ 

j(0p,0p)GWp 



n (i/0 d 



> 



E w ({ J p» 



{jp> 



(A) 

A_ 

4^j 



2\ 



x exp [-pE* ({J p }) - /3L d A max A] 



> max exp 



/3£*({J p })+lnW({J p }) 



L\ d . A 2 



In — - /3£ d A max A 

47T 



(41) 



and obtain 



Lim— t 

L d 



H C \{J P + 1, 9 P , </> P }) - H C \{J P , 9 P , P }) 
2sl d + 1 



< Lim- 



2s 2 l 2d 



A max = 0. (35) 



Since A > is arbitrary, we choose A = I d . We then 
find 



Lim/ < Lim inf —? 

{S P } L d 



1 



H cl ({J p ,9 p ,<j> p }) - -lnW({J p }) 



(42) 
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From the upper bound (42) and the lower bound (36), 
we obtain the desired variational expression for the free 
energy, 



Lim/ 



Lim inf — r 



1 



H c \{J p ,d p ,cf> p })--lnW({J p }) 



(43) 

The reason why we can reach the variational expres- 
sion (43) is that the degrees of freedom to be considered 
is dramatically decreased by the coarse graining. 

We can express the free energy by the continuous func- 
tional. If we see the system in the length scale so the size 
of the system is unity, 



inf 

5(o=) 



d d x / d d yU(x-y) 
[os) d J[o,i) d 



]T X a S a (x)S a (y) 



a—x,y,z 



-h- [ S{x)d d x- \ [ a{S{x))d d x 
J[os) d P J[as) d 

= miT[S(x)}. (44) 

Here we defined the scaled potential 

U(x) = Um(-fL) d <j)(-fLx), 
and the entropy function 

a(S(x))=Um^\nW(J p ), 

with S(x) = S p /(sl d ) and x = Lim(r p /L), where r p is 
the central position of the cell C p . Notice that the length 
of the spin vector does not exceed unity, [S^a;)! < 1. 

The free energy functional ^[^(x)] consists of two 
parts, 



F[S(x)} = U[S{x)] - - [ a(S(x 
P J[o,i] d 



))d d x, 



where 



U[§\ = -\( d d x f d d yU(x-y) 
z J[o,i) d J[o,i) d 

V \ a S a (x)S a (y)-h- ( S{x)d d x 

a=x,y,z • / [0,l) d 



(45) 



is the internal energy functional. It corresponds to the 
energy of a given configuration S(x)). j cr(S(x)d d x cor- 
responds to the entropy of a given configuration S(x). 
Because many microscopic states correspond to a sin- 
gle coarse-grained configuration S(x), this entropic term 
arises. Thus, the calculation of the free energy is re- 
duced to the minimization problem of the free energy 
functional. 



It is noted that if the uniform spin configuration, 
S(x) — m independent of x, is assumed, the free en- 
ergy is exactly equal to that of the mean-field model (the 
model with 4>(x) = 1/N) and independent of the poten- 
tial (j)(x); 



J"[5(a;) = rfi] = /mf(/3, m, h). 



(46) 



The explicit form of fMF(P,m,h) is given in Eq. (47). 
This observation implies that we can replace the interac- 
tion potential 4>{x) by the simple infinite range potential 
1/N as long as the magnetization profile is uniform in 
equilibrium. While, the spin configuration must be inho- 
mogeneous whenever this replacement is not allowed. 

Finally, I make a few remarks on the relation to pre- 
vious works. Gates and Penrose [20, 21] obtained the 
similar variational expression of the free energy. How- 
ever, the free energy functional obtained in the present 
paper [Eq. (44)] differs from that in those works with re- 
spect to the length scale of the free energy functional. 
In the present work, the position variable x runs over d- 
dimensional unit cube, x € [0, l] d , and the system size is 
set to be unity. It enables us to treat the long-range inter- 
action potential like l/r a , a < d. On the other hand, in 
the free energy functional obtained by Ref. [20] , the posi- 
tion variable x runs over the whole d-dimcnsional space, 
x e M. d . In that case, we cannot treat the long-range 
interaction like l/r a with a < d because the integrated 
value of the potential is divergent. If we consider only the 
van der Waals limit, however, the free energy functional 
in Ref. [20] has an advantage that it has information on 
the mesoscopic length scale. In the expression (44), we 
cannot obtain any information on the length scale ~ 7 , 
e.g. the form of the interface in an ordered state with 
phase separation. 

In the work by Kiessling and Percus [22], the concep- 
tually different but mathematically equivalent limiting 
procedure was introduced to study the liquid-vapor in- 
terface in a many particle classical system. That is, they 
took the limit of infinitely many particles in a finite do- 
main. As a result, they derived the similar variational 
expression as Eq. (44) on the thermodynamic potential. 
It shows the similarity between bulk properties of a sys- 
tem with a nonadditive long-range interaction and lo- 
cal properties of a system with a finite-range interaction 
and with a very large number of particles in a bounded 
domain. In the works of Refs. [10-12], from the vari- 
ational expression (44), simple sufficient and necessary 
conditions for the exactness of the mean-field theory are 
derived in classical spin systems, which is presented in 
the next section. 



VI. EXACTNESS OF THE MEAN-FIELD 
THEORY 

Since the calculation of the free energy is reduced to the 
minimization of the free energy functional expressed by 
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the classical variables, we can apply the same argument 
as the case of classical spin systems. In this section, we 
briefly review the result of the variational expression (43) 
or equivalently, the continuum expression (44). 

Hereafter, we assume periodic boundary conditions 
and apply the nearest image convention to pair inter- 
actions. 



A. Canonical ensemble 

We assume that the system is in contact with a ther- 
mal reservoir at a temperature T = and there is no 
conserved quantity. In this case, the free energy is inde- 
pendent of the interaction potential 4>(-) and it is equal 
to the mean-field free energy 



/mf(/3) 



a— x,y,z 

= mmf M F(P,m,h) 



\ x aml - h ■ m - ^a(m) 



(47) 



where 



m 



( S(x)d d x. 
J[o,i) d 



This result is called "exactness of the mean-field theory" 
in long-range interacting systems. 

It is convenient to consider the canonical ensemble with 
the fixed magnetization, although the magnetization m 
is not a conserved quantity. Roughly speaking, the free 
energy 

f(/3, m, h) = min F[S(x)] (48) 

{S(x):m=f [0A)d S(x)d d x} 

is related to the probability P{rh) of the states with 
the magnetization rfi by P(rfi) <~ cxp[— /?/(/?, m, h)}, al- 
though the meaning of "the states with the magnetization 
m" is ambiguous in quantum spin systems because of the 
quantum fluctuation. This function is well-defined and 
indeed useful when we analyze the case where there is 
some conserved quantity like the energy in Sec. VI B or 
m z in the XXZ Heisenberg model discussed in Sec. VII. 

Later, we write f((3,rn,h) simply as /(/3,m) because 
the dependence on h is trivial. 

In the van der Waals limit, the scaled potential is the 
S function, U(x) = S(x), and in this case we can show 



fW,m) = f^ F (p,m) 



This corresponds to the Lebowitz-Penrose theorem. 

In the nonadditive limit, U{x) = (f>(x). The Fourier 
expansion is given by 



U(x) = U ™ e 



2-Kin-x 



(50) 



(51) 



If we define 

CAnax = max U n , 

n£Z d \0 

we obtain [10, 11] 

/MF(/3, m) ~ U maK Af M F{/3U max ,m) 

</(/?, m) </M F (/3,m), (52) 
where A/ M f = /mf - /mf. and moreover 

/(P)</mf(P) (53) 

when the matrix 



m)] ab 



d 2 h 



MF 



dm a dmb 



WU max ,m) (54) 



is non-positive. 

From the inequality (52), we find that 
/(/3,m) = /mf(/3,to) in the parameter (/3, m) such 

that fMF{PU ma x,m) = /M F (^max,m). 

We define by 

P = sup [p > : f MF (i3,m) = fuFiP,™) Vtti] . (55) 

It is then concluded that /(/3,m) = /mf(/3, m) for any 
m when /3J7 max < /3. This plays important roles in the 
following sections. 



B. Microcanonical ensemble 

In this section we consider the microcanonical ensem- 
ble without any other conserved quantities. The case 
in which there is another conserved quantity like m z is 
treated in Sec. VII. The microcanonical entropy is ob- 
tained by 



sup 

{£(*)} 



J[o,i) d 



u{S{x))d d x : e = U[S(x)] 



, (56) 



where the internal energy functional U[S(x)] is defined by 
Eq. (45). This expression is common to classical spin sys- 
tems, so the same argument as classical spin systems [12] 
can be applied. 

We can easily verify 



smf(e) < s(e) < %f( £ )i 



(57) 



where smf is the microcanonical entropy in the mean- 
(49) field model, <j>(x) = 1, and Sy F (e) is the concave enve- 
lope of smf(e) or the function obtained by performing the 
Legendre transformation twice on Smf(e)- From the in- 
equality (57), we can conclude that the mean-field theory 
is exact in the microcanonical ensemble, s(e) = Smf(e), 
if the canonical and the microcanonical ensembles are 
equivalent in the mean-field model because smf = sJf F 
in this case. 
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In the case where the canonical and the microcanonical 
ensembles are inequivalcnt in the mean-field model, more 
sophisticated treatment is necessary. In this case, it was 
shown in [12] that 



smf(e) < s(e) 



< min < SM F (e), sup 

I S'>6 



smf(e') + 



-(e-e') 



(58) 



Here /3 is given by Eq. (55). Moreover, it was shown that 

a(e) > s MF (e) (59) 

when the matrix £(/3MF(e)f/max, rh e q(e)) is nonpositive. 
The definition of the matrix C is given by Eq. (54). The 
quantity /3mf(e) = dsyiF/de is the inverse temperature 
of the mean-field model at the energy e, and rh eq (e) is 
the equilibrium magnetization of the mean-field model, 



m eq (e) ee argsup 



a(A):e = ~ ^m 2 a - h ■ 



in 



(60) 

The result of the inequality (58) is as follows. We con- 
sider the case that smf(e) < SM* F (e) in the region of 
e a < e < Sb- We define 



= sup 



dfMg(f) 
de 



: e e (e a ,£ b ) 



(61) 



If /3*J7 max < P, then s(e) — smf(e) holds for all the val- 
ues of e. Therefore, the states with the negative specific 
heats in the mean-field model are maintained for sys- 
tems with general long-range interaction <p(-) as long as 
the condition f3*U max < f3 is satisfied. 

On the other hand, the inequality (59) means that 
when the matrix £(/3MF(e)^max, m eq (e)) is non-positive, 
the equilibrium state in the mean-field model, S(x) = 
rh eq (e), becomes unstable and the inhomogencity ap- 
pears. In this case, the mean-field theory is not exact. 



VII. EXAMPLE: THE SPIN-1/2 XXZ MODEL 

As an application of the theoretical framework, we 
consider the spin-1/2 XXZ model, whose Hamiltonian 
is given by Eq. (1) with X x = X y = X± > and X z > 0. 
The Hamiltonian is given by 

H = -\H W(7r«) [A±(*f s? + s y s y) + x zS f S |] . 

ij 

(62) 

The mean-field counterpart of this model [7*^(7^) — 
1/N in Eq. (62)] was studied by M. Kastner [4, 5]. 

When the magnetization is not fixed, the exactness 
of the mean-field theory always holds in the canonical 
ensemble, f((3,h) = /mf(/?> which is a general re- 
sult from Sec. VIA. Moreover, in this model, the mean- 
field theory is exact even in the microcanonical ensemble, 



s(e) = smf(e)- It is resulted from the fact that the micro- 
canonical entropy of the mean-field model corresponding 
to Eq. (62) is a concave function of e. It implies that 
Smf(e) = Smf( £ ) an d> thus, s(e) = smf(e) from Eq. (57). 

If the system is isolated, the z component of the mag- 
netization m z is conserved. In this section, we consider 
the canonical and the microcanonical ensembles with a 
fixed value of the z component of the magnetization m z . 
In this case, as we will see below, the result is nontrivial 
and the exactness of the mean-field theory is violated for 
certain parameter region. 



A. Canonical analysis 

First, we investigate the model in the canonical ensem- 
ble. The free energy as a function of (3 and m z is given 

by 

fifi, m z )= inf /(/?, m) = f(J3, m eq (/?, m z )). (63) 

{m x ,m y } 

Of course, {rh cq (/3,m z )) z = m z . 

The analysis of the previous sections leads us to the 
following result. From the inequality (52), we obtain 

/ MF (^,TOc q (^,TO^)) - U max Af M F(PU mm ,m e q(/3,rn z )) 
< f(P,m z ) < f M F(/3,m z ). (64) 

Therefore, if /??7 max < /3, f(P,m z ) = fMF{P,m z ) holds 
because AfMF{flU max , rh) = 0. Moreover, if the ma- 
trix C(/3U majX , m eq (/3, m z )) has negative eigenvalues, the 
mean-field model is not exact, f(/3,m z ) < /mf(/?, Tn z ). 

Since the entropy function of the spin-1/2 system is 
given by Eq. (A4) in Appendix A, the mean-field free 
energy with a fixed m is given by 

/mf(/3, m) = ~ [A_l(to^ + my) + X z m 2 z ] 

i/i + N i + H + i_H ln i_^y (65) 

V 2 2 2 2 J V ' 

In order to obtain the mean-field free energy as a function 
of m z , it is necessary to minimize /mf(/3,w) for m x and 
m y . By doing that, we find in the mean- field model 

(jneq)x = (iTL e q)y = when m z > tanh(/3Aj_m z ), 
|m eq | = tanh(/3Aj_|m eq |) when m z < tanh(/3Aj_m z ). 

By differentiating twice the mean-field free energy 
/mf(/3,w), we obtain the Hessian matrix C(f3,m) (see 
Eq. (54)). We can obtain /3, which plays important roles 
to judge whether the mean-field theory is exact at a 
given parameter (j3,m z ), by analyzing the eigenvalues 
of £(/3, rh). Actually, we obtain 

= min {AI 1 , A; 1 }. (67) 

See Appendix C for the proof of Eq. (67) . 
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By the reason mentioned in later analysis on the mi- 
crocanonical ensemble, f((3,m z ) = fMF{fl, m z) when 
Aj_ > A z . The result is non-trivial when A z > Aj_. In 
this case, = 1/X Z , and the mean- field model gives the 
exact free energy, f((3,m z ) = /mf(/3, ?^ z ) at least for 

P < l/(t/maxA z ). 

By calculating the minimum eigenvalue of 
£(/3C/ max ,m eq (/3)), we obtain the sufficient condi- 
tion under which the exactness of the mean-field theory 
does not hold. 

In Fig. 1, we plot the two regions in (T = m z ) 
plane. One is the region which is given by j3 < 
l/(U max X z ) (black region). In this region, the mean-field 
theory is exact. The other is the region where the ex- 
actness of the mean-field theory is violated (meshed gray 
region), which is given by the condition that the min- 
imum eigenvalue of £((3U max , m eq (/3)) is negative. The 
parameter ?7 max is set to be 0.3. This value corresponds 
to the case that <fi(r) oc 1/r in the two-dimensional lat- 
tice. In the white region in Fig. 1, we cannot determine 
whether the mean-field theory is exact from our analy- 
sis. However, in the white region, the equilibrium state 
predicted by the mean-field model is at least locally sta- 
ble because all the eigenvalues of £((3U maK , m eq (/3)) are 
positive. 



B. Microcanonical analysis 

Next, we consider the model in the microcanonical 
ensemble. The mean-field entropy was obtained and 
the property of the mean-field model was studied ex- 
tensively [4, 5]. Some of the results obtained here will 
overlap these studies. The microcanonical entropy in the 
mean-field model is given by 



When m z is conserved, the analysis given in Sec. VI B 
must be slightly modified. Inequality (58) is altered to 



smf(e, m z ) 
where 



1 + 9 , 1 + 9 1 - 9 , 1-5 f(tQ , 
In— —In——, (68) 



g = y/m>{l - X z /X ± ) - 2e/X±. 



(69) 



Because there is the conserved quantity m z besides the 
energy e, the domain of the entropy is not necessarily a 
convex set. The domain is given by 



s M F(e,m z ) < s(e,m z ) 
< min < s^ F (e, m z ), sup 

I £'>£ 



SMF{e',m z ) + — (e - e') 



(71) 



where s^ F (e,m z ) is the function obtained by applying 
twice the Legendre transformation to SMF(E,m 2 ) with 
respect to both variables (e,m z )[23\. From the inequal- 
ity (71), it is concluded that s(e,m z ) — Smf^,^) at 
least when (e,m z ) E T> M f and /?mf(£, m z ) < $/U ma , x . 
Here, Pmf{s, rn z ) = ds M F{e,m z ) / de. 

The condition of Eq. (59) is not altered; if the matrix 

C(Pmf(£, m z )U maK , m oq (e, m z )) 

has a negative eigenvalue, it is concluded that s(e, m z ) < 
smf(e, m z ). 

When A^ > X z , the mean-field entropy smf(£, rn z ) is 
a concave function of e and m z , see Eq. (68). Therefore, 
because of the inequality (71), smf(£,w z ) = s^ F (e,m z ) 
in this case. Since the canonical free energy is obtained 
by the microcanonical entropy, it is also concluded that 
f(f3,m z ) = /mf(/3,to z ) when Aj^ > X z . The exactness 
of the mean-field theory holds when Aj^ > A z both in 
the canonical and in the microcanonical ensemble with a 
fixed value of the magnetization m z . 

On the other hand, the result is non-trivial when A z > 
Aj^. In Fig. 2, we plot the region of Vmf (inside the black 
solid line), /3 M F{£,m z ) < /3/C/ max = l/(U max X z ) (black 
region), and the region where 

£(I3mf(£, m z )C/ max , m cq (e, m z )) 

has a negative eigenvalue (meshed gray region), in the 
(e,m z ) plane. In the black region, the mean-field the- 
ory is exact. In the meshed gray region, the mean-field 
theory is not exact. In the white region inside the black 
solid line, we cannot determine whether the mean-field 
theory is exact from our analysis, although the equilib- 
rium states predicted by the mean-field model is at least 
locally stable. 



V 



MF 



(e,m z ) : £ < — -m z 

a ^ Aj_ - A z 2 1 
and e > mi - -X± 



(70) 



The domain of the model (62), which is referred to as V, 
and that of the mean-field model, £>mf, differ in general, 

V D Dmf- 

Exactness of the mean-field theory is trivially violated 
for parameters (e, m z ) G 2?\2?mf- 



VIII. SUMMARY AND DISCUSSION 

In this paper, we analyzed quantum spin systems with 
long-range interactions. We reduced the calculation of 
the free energy to the minimization of the free energy 
functional, which is expressed only by the classical vari- 
ables. We obtained some results from this variational 
expression of the free energy. One of them is that the 
free energy is always identical to that of the mean-field 
model when there is no conserved quantity. When there 
is a conserved quantity such as m z , the mean-field model 
is exact at least for /3i7 max < j3. The key quantity J7 max 
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(a) (b) (c) 




0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 

T T T 

FIG. 1: The abscissa is the temperature T = 1//3 and the ordinate is m z . In the black region, the mean-field theory is exact, 
and in the meshed gray region, the mean-field theory is not exact. We cannot determine whether the mean-field theory is exact 
in the white region. The parameters A z and I7 max are set to 1 and 0.3, respectively, (a) Ax = 0.2, (b) Ax = 0.28, and (c) 
A ± = 0.4. 




FIG. 2: The abscissa is the energy per spin e and the ordinate is m z . The region of 2?mf is inside the black solid line. 
In the black region, Pmf(s, m z ) < /3/t/ m ax = l/(I/maxA z ) and the mean-field theory is exact. In the meshed gray region, 
'C( / 0MF(e,m 2 r)?7 m ax,meq(e, m„)) has a negative eigenvalue, and the exactness of the mean-field theory is violated. We cannot 
determine whether the mean-field theory is exact in the white region inside the solid line. The parameters A z and Umax are set 
to be 1 and 0.3, respectively, (a) Ax = 0.2, (b) Ax = 0.28, (c) Ax = 0.4, and (d) Ax = 0.6. 



is determined by the form of the interaction potential 
<f>(jr) and the spatial dimension d. Moreover, if the ma- 
trix £(/3£/ max , m eq (/3, m z )) has a negative eigenvalue, the 
state with the uniform magnetization profile predicted 
by the mean-field model is locally unstable, and inhomo- 
geneity which is not predicted by the mean-field model 
arises. These results are the same as those in classical 
spin systems with long-range interactions. 

We analyzed the systems also in the microcanonical 
ensemble, and obtained the same results as those in clas- 
sical spin systems. When the magnetization is not con- 
served, the microcanonical entropy is identical to that of 



the mean-field model as long as the canonical and the 
microcanonical ensembles are equivalent in the mean- 
field model. However, if these two ensembles are in- 
equivalent in the mean-field model, the microcanonical 
entropy may depend on the interaction potential 4>(x). 
The obtained results in this case are as follows. The 
mean- field model is exact if /3MF(£)U max < $ for all e 
such that smf(e) < Smf( £ )- On the other hand, when 
£(,0MF(£)kmax, rn cq (e)) has a negative eigenvalue, the 
microcanonical entropy is not identical to that of the 
mean-field model because the uniform state predicted by 
the mean-field model is locally unstable. The analogous 
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result holds when the magnetization m z is conserved. 

As a concrete example, we investigated the spin-1/2 
XXZ model with a long-range interaction. When the 
magnetization is not conserved, the result is trivial (the 
mean-field theory is exact) both for the canonical and 
the microcanonical ensembles. Therefore, we studied the 
model with a conserved m z . As a result, it is shown that 
the exactness of the mean-field theory is violated in a 
certain parameter region for both the canonical and the 
microcanonical ensembles. 

In cold atom experiments, we can set up the system 
which is well isolated from the environment, but the z- 
componcnt of the magnetization is not a conserved quan- 
tity in general. Such a system will be described by the 
microcanonical ensemble without restriction of the value 
of m z . In this ensemble, the spin-1/2 XXZ model with 
long-range interactions is equivalent to the correspond- 
ing mean-field model, as mentioned above. If we seek to 
observe the violation of exactness of the mean-field the- 
ory in long-range interacting systems, it is necessary to 
consider another model. 

In order to observe the violation of the exactness of 
the mean-field theory in the microcanonical ensemble, 
the mean-field model corresponding to the system must 
satisfy smf(e) 7*= s mf( £ )- This condition means that the 
canonical and the microcanonical ensembles are inequiva- 
lent in the mean-field model [24] . Therefore, if the actual 
interaction is not infinite-range (or Curie- Weiss type), 
not only the ensemble inequivalence predicted by anal- 
ysis of the mean-field model, but also the violation of 
the exactness of the mean-field theory will be observed. 
Although it is known that the mean-field models repro- 
duce many properties of general long-range interacting 
systems even quantitatively, the infinite-range (or Curie- 
Weiss type) interaction does not become an adequate 
idealization of actual long-range interactions in such a 
situation. 

Kastner pointed out that the inequivalence of the 
canonical and the microcanonical ensembles can be ob- 
served in a long-range interacting system undergoing a 
temperature-driven first order phase transition [4]. A 
typical example of such spin models which possibly ex- 
hibits the violation of the exactness of the mean-field 
theory is the spin-1 Ising model with anisotropy, 



(72) 



The corresponding mean-field model, which is ob- 
tained by replacing j d <j)(jrij) by 1/iV, undergoes the 
temperature-driven first order phase transition for cer- 
tain range of the parameter D. However, this model is a 
special one in the sense that it is a classical spin system. 
It is expected that more general quantum spin systems 
with the anisotropic term, 



behaves similarly, although it is not proven straightfor- 
wardly from the formalism given in this paper. The pres- 
ence of the nonlinear external field f Xii( s i) 2 makes it 
difficult to treat by the formalism given in this paper. A 
future work should treat nonlinear external fields such as 
£> J7(sf) 2 , which will promote an experimental realiza- 
tion of characteristic properties of long-range interacting 
systems. 
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Appendix A: Calculation of W(J P ) 

We seek to calculate the weight function W(J P ). From 
the definition, W(J P ) is the total number of trajectories 
with length l d from to J p . We can obtain the expression 
of W(J P ) by counting such trajectories, but we can do it 
more easily. After all, W(J P ) is the number of states with 
the total spin J p and the total magnetization M p = J p . 



The states with M p — X have J p 
States with J p > X 



X or J p > X 



1 and M p — X are obtained by 



multiplying S~ to states with M p = X + 1. Therefore, if 
we refer to Wm{X) as the number of states with M p = X, 
the weight function is given by 

W{X) = Wm(X)-W m {X + 1). (Al) 

For instance, in the case of spin-1/2, Wm(X) is given 



by 



W M {X) 



l d \ 



(*+x)\{$-x)\ 



(A2) 



Therefore, 



W{J P ) = W M (J P ) - W M (J P + 1) 
2J V + 1 l d \ 



Applying the Stirling formula, we obtain 



(A3) 



a(S) = Lim^ In W(J P ) 

_i±M m i±M_i_M ln i_M (A4) 
2 2 2 2 v ; 

The weight function and the entropy function are cal- 
culated similarly for higher spins (,s = 1,3/2,...), though 
the calculation becomes complex. In the case of s = 1, 
the result is 



H = — 



—y 

2s 2 ^ 



i a -x 



I" 



a—x,y,z 



(73) 



N=0 



[(l d -N + X)/2}l[(l d - N - X)/2]\N\ ' 



(A5) 
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and 



a(S) 



where 



2 2 

1 -x- \S\ 1 -x- \S\ 



— x\nx, (A6) 



-1 + 



Appendix C: Proof of Eq. (67) 

In the model described by Eq. (62), the mean- field free 
energy is given by Eq. (65). The Hessian matrix £(/3, rh) 
defined by Eq. (54) determines the convexity of the mean- 
field free energy. If the minimum eigenvalue of £(/?, rh) 
is positive for all rh, the mean-field free energy is convex 
for rh and /mf(/3, rh) — /m*f(/3>to). Therefore, j3 corre- 
sponds to (i such that the minimum eigenvalue of £(/3, rh) 
becomes zero at some rh. 

Eigenvalues of £(/?, rh) are calculated and given by 



Appendix B: Proof of Eq. (40) 

We derive the inequality (40). For simplification, we 
write e £ u p (A) if we can express the unit vector e by 
e = (sin 9 cos 4>, sin 9 sin 0, cos 9) with (0, $) e u p (A). 

First, we show that \a — b\ < 2 A for arbitrary unit 
vectors a, b £ co p (A). We have 

\a-b\ = \{a-e* p )-(b-e* p )\ < \a-e* p \ + \b-e* p \. (Bl) 

We choose the z axis as the direction of the vec- 
tor e* p . Any x £ uj p (A) is represented as x — 
(sin 9 cos (j>, sin 9 sin 4>, cos 9) with < 9 < A, which yields 
\x-e* p \ = 2sin(0/2) < 9 < A. Therefore, from Eq. (Bl), 
we have 

\a- b\ < \a- e* p \ + \b - e* p \ < 2A. (B2) 
We now show Eq. (40). From the definition, we have 

sup H cl ({J p ,6 p ,<l> p })-E*({J p }) 

(e p ,<p p )eu p (A) 



sup 

(9 p ,<p p )eu p (A) 



! (L/l) d 
7^2 1 d< t>P1 J P J <} 



U(e P ) a (e q ) a (e* P ) a (e%) a ) 



(L/iy 



< ^ SUp ^ ^ S lUpqJpJq^ 



(e p ,0 p )GWp(A) 



P,9 



(B3) 



Here we define e p = (sin # p cos </> p , sin p sin (/> p , cos 6* p ) and 
e* p = (sin 9* cos 0*, sin 8* sin 0*, cos (j) p ). Because J p < 
sl d and 



|6p • 6^ 



< |e p - e* p \ + |e 9 - e* 

< 2A, 



(B4) 



we obtain the inequality (40), 

sup H c \{J p ,9 p ,0 p })-E*({J p }) 

(e p ,4> p )euj p (A) 

< A max Al 2d ]T 7 % 9 = A max AL d . (B5) 



i fa ^ \ 1 1 i 1 + H 

/±(/? '™ ) = 2^R ln r 



(Cl) 



to 



TO 



+ 



1 



— Aj_ -A z + B|to| 



± V(A± - A z + B|to| 2 ) 2 + 4(A Z - X^)B(\rh\ 2 - to 2 ) , 

(C2) 

where 



B EE 



1 1 

2/3 "H 



1 , l + \m\ 
In ■ 



1 — TO TO 1 — TO 



By using these eigenvalues, /3 is given by 



= sup | 



p > : min 



inf Zo(/3 ; to), inf l-(/3, rh) 



Obviously, we find 



inf l (/3,rh) = + -. 

m p 



>0|. 

(C3) 

(C4) 



When A^ > A z , it is shown that rh) > lo{(3,rh). 

Therefore, f3 = 1/Aj_ in this case. 

Next we consider the case of X z > X±. We rewrite 
rh) as 

/ fa ^ 1 1 . 1 + |to| 



1 

+ 2 



2/3\rh\ l-\rh\ 

\ ± -X z + B\rh\ 2 



- V(*±- - *z - B\m\ 2 ) 2 - 4(A, - \x_)Bml . 

(C5) 

From this expression, 



^m)>-A 2 + -L^lnl±^[ 

2/3 TO 1 — TO 



Hence, we have 



inf l_ (p, m) >-\ z + -. 

m p 



(C6) 



(C7) 
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On the other hand, from Eq. (C2), From Eqs. (C7) and (C9), we obtain 

l-W,m) < ^^ln 1 + |"!| -X z + B\m\ 2 i 

V ; -2/3|m| l-|m| 11 inf Z_(/3,m) = —\ z + \. (CIO) 

11 m p 

= -A» + 7 n FHa - (C8) 
pi — |m| J 

Thus we have Therefore, when A z > Aj_, we have /3 = 1/A Z . 

The analysis in this appendix concludes that (3 = 
Ml_(f3,m)<-\ z + -. (C9) min(l/A ± ,l/A z ). 

m p 
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